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Abstract 

We analyse the combustion process in an spark ignition engine using the experi- 
mental data of an internal pressure during the combustion process and show that 
the system can be driven to chaotic behaviour. Our conclusion is based on the 
observation of unperiodicity in the time series, suitable stroboscopic maps and a 
complex structure of a reconstructed strange attractor. This analysis can explain 
that in some circumstances the level of noise in spark ignition engines increases 
considerably due to nonlinear dynamics of a combustion process. 
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1 Introduction 



It is known that the cyclic combustion variability is one of the main charac- 
teristics for spark ignition (SI) engines. If cyclic variability were eliminated, 
there would be even 10% increase in the power output of the engine [1]. Cyclic 
variations in the combustion process of SI engines have been a subject of an 
intensive research in last 40 years. Heywood [2] identified three main factors 
influencing cycle-to-cycle variations: aerodynamic in the cylinder during com- 
bustion, the amount of fuel, air and recycled exhaust gases supplied to the 
cylinder and a mixture composition near the spark plug. 

Recently, in the context of engine control, there appeared a lot of papers 
on dynamic phenomenon identifications and predictions of engine behaviour 
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[3,4,5,6,7]. Many of them discuss the problem of a high noise level which make 
the engine difficult to control. On this basis several stochastic models have 
been applied to pressure cyclic variation analyzing [4,5]. 

As the combustion of gases in the cylinder may be viewed as a nonlinear 
dynamical process [1,3], it may be studied by tools belonging to the theory 
of nonlinear systems. The variations might originate in nonlinear dependence 
of the peak cycle temperature and pressures on the initial conditions at the 
beginning of compression due to exhaust recirculation and a mixture prepa- 
ration process. Such nonlinear approach has been already started by Daw et 
al. [3] where he suggested the possible chaotic nature of combustion. Our pa- 
per explanation is going in the same direction, however in our approach we 
adopted the novel measurement method basing on optical fibers to measure 
the pressure directly [8] . 



2 Chaotic oscillations of internal pressure 

We start our analysis monitoring the measured time series for various spark 
advanced angles Aa^. The crankshaft frequency was chosen to be © = 14.17 
Hz (© = 850 RPM) as an idle speed of the engine. In a combustion process 
only one of each two cycles of crankshaft rotation are coincided with com- 
bustion leading to pressure oscillation with a period 2/0 = 0.141s. Note that 
depending on ignition timing (Actx) the system behaves differently. Starting 
with AcKj = 5° (Fig. la) we observe initial periodic oscillation. Increasing 
Aq;^ to 20° the situation docs not change very much nevertheless, one can 
note that in the background of a periodic motion some singular instabilities 
occur (i.e. at t=2.8s at Fig. lb) after a long interval of a stable motion (Fig. 
lb). Interestingly, going further to a larger value of the spark advanced angle 
Aoij; = 30° changes the system behaviour completely. The corresponding time 
history is plotted in Fig. Ic. In this case the system has lost its initial periodic 
nature undergoing to a more complicated type of motion. In the last panel 
(Fig. Id) we show the power spectrum for this interesting case. The continues 
spectrum indicates that system can be in a chaotic state. 

To clarify our conjecture about chaotic combustion we decided to present re- 
sults on the stroboscopic map. In this aim we prepared Fig. 2 where we plotted 
the internal pressure coincided with the position of the crankshaft rotation an- 
gle. For enough log time history we see a kind of collection of Poincare maps 
plotted simultaneously for large number of performed measurements. Namely, 
our crankshaft was divided into 512 angles. These gives 1024 measurement 
points per two cycles of a crankshaft rotation corresponding to one cycle of 
combustion. One can see that initially seen periodic motion (Fig. la) has 
already some modulation of combustion slope. The critical region of modula- 
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tion was however located behind the maximum pressure peak leaving its value 
p„„,,.(<6) practically unaffected. Obviously, this modulation, rather difficult to 
observe in Fig. la, introduces another time scale with longer period. Note that 
for larger value of Aa^ (Aa^ = 20") that critical region has moved towards 
the maximum pressure peak giving rise to observable changes of the maxi- 
mum pressure value (Fig. 2a). These changes become the most dramatic for 
= 30" where the the modulations of pressure exceeds its average value 
(Fig. 2c). To illustrate the final effect of pressure instabilities we show, in Fig. 
2d maxima of pressure Pmax in sequential cycles n. In this figure squares corre- 
spond to — 20" and circles to 30°, respectively. Note that the maximum 
value of pressure Pmax for Aq;^ = 5° would be a constant line in the scale of 
Fig. 2d. 

To this end we plot in Fig. 3 the reconstruction of the chaotic attr actor. 
Changing the parameter Aa^ which can be regarded as a bifurcation parame- 
ter we see that the attractor which was in principle two dimensional enlarge its 
dimensionality with increasing Aa^- The expected minimal embedding dimen- 
sion was chosen arbitrary to be equal to 3 as on the basis that it is the smallest 
dimensionality enabling chaotic solutions. The characteristic time delay value 
T = 0.138s was chosen simply as a number smaller than the characteristic 
period in the system 2/0 = 0.141s (Fig la-c). Note 2/0 as a period of a 
parameteric excitation is uneffected by changes of the bifurcation parameter 
Aaz- Starting with basically 2-d structure of the attractor in tree dimensions 
(Fig. 3a) one can follow small changes of it for Aaz — 20" (in Fig. 3b) and a 
large quahtative change in Fig. 3c for Aa^ — 30°. 



3 Conclusions 

In this paper we examined the oscillations of combustion internal pressure. It 
appeared that, in some conditions, intermittency is capable to drive the system 
into chaotic region. We have shown that the change of a spark advance angle 
Aa^ is making a significant effect on the combustion dynamics. It is clear that 
it influences directly the flame initiation phase which occurs on a larger time- 
rate for an increasing spark advance. In this way the cyclic variations become 
an inherent phenomenon in time and in space. Chaotic nature of combustion 
in SI engines appears in strong sensitivity on such conditions. It is worth to 
note that Figs. Ic and Id, for Aa^ = 30° can be easily associated with a 
chaotic process, but Figs, lb and 2b show that even for the advance angle 
Aaz — 20° main instabihties appear. This fact leads us to conclusion that 
the mechanism of transition to chaos should be based on the intermittency 
phenomenon [9,10,11]. 

Our preliminary results indicate that the noise SI engines ought to be con- 
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nected with nonlinear dynamics. Moreover analysing the combustion process 
for different system parameters we observed similar behaviour. I should be 
noted, however that in parametric systems with an additional self-excitation 
the problem of distinguishing between chaotic and multi frequency quasi- 
periodic motions is not an easy task [12]. To tell more about dynamics of 
that particular combustion method one have to use more sophisticated non- 
linear methods [13]. We plan to conduct such analysis in a future report [14]. 
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Figure Captions 



Fig. 1. Time series of pressure in a combustion process the crankshaft frequency 
= 850 RPM (a-c) for various spark advance angles Aa^. The power spectrum for 
Aa^=30° (d). 

Fig. 2. Stroboscopic projection of internal pressure p{t) on a relative position of a 
crankshaft in a combustion process for various Aa^ (a-c). Maxima of pressure 
Pmax in sequential pressure cycles n where squares correspond to Aa^ = 20" and 
circles to 30", respectively (d) (squares are connected by a dotted line for better 
clarity) . 

Fig. 3. Reconstruction of the attractor for a combustion process in a SI engine. 
Characteristic time delay was assumed to be r = 0.138 s. 
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